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We demonstrate that the system of locahzed magnetic mo- 
ments in HoNi2B2C can be described by the 4-positional clock 
model. This model, at a proper choice of the coupling con- 
stants, yields several metamagnetic phases in magnetic field 
at zero temperature in full agreement with the experimental 
phase diagram. The model incorporates couplings between 
not nearest neighbors in the direction perpendicular to the fer- 
romagnetic planes. The same model leads to a c-modulated 
magnetic phase near the Curie temperature. The theoret- 
ical value of the modulation wave-vector agrees surprisingly 
well with that observed by the neutron diffraction experiment 
without new adjustable parameters. 



In the works by Rathnayaka et al. ||l| and by Canfield et 
al. Q transport and magnetic measurements on HoNi2B2C 
for various magnetic fields and low temperatures have been 
reported. The magnetic phase diagram for HoNi2B2C with 
fields in the a-b plane is of particular interest. 

In this compound easy magnetization axes are identified 
with crystallographic directions (110) and (110). The low 
temperature magnetization data show the existence of 4 meta- 
magnetic phases. The low field phase has been identified by 
neutron difi'raction experiments and magnetic measure- 
ments Q with the antiferromagnetic phase, which we denote 
symbolically ][. The phase boundaries and magnetization in 
other phases versus magnetic field found in the experiment |^] 
can be readily explained by assuming that the remaining three 
phases are as follows: phase 2 - ]][, phase 3 - TT— and the 
high-field phase 4 - t. It means that | of the spins in the 
phases 2 and 3 are parallel to one of the easy axes whereas the 
remaining i is antiparallel and perpendicular, respectively, to 
the same axis. Note that all metamagnetic phases are stoi- 
chiometric in the meaning that the concentrations of spins 
parallel, antiparallel, or perpendicular to the reference axis 
are rational numbers. The phase diagram of HoNi2B2C at 
zero temperature is especially simple if the components of 
magnetic field H^, Hy are chosen as variables. The experi- 
mental phase diagram of HoNi2B2C is shown in Fig. 1 (in the 
original work it has been presented in polar coordinates h, 
Oh). 

Siegrist et al. and Huang et al. ^ determined the struc- 
ture of Lu and Ho 1:2:2:1 compounds as the body-centered 
tetragonal lattice with the space group M/mmm. The x- 
ray structure analysis and the neutron-scattering experiments 
performed by Goldman et al. and Grigereit et al. 
showed that an incommensurate modulated magnetic struc- 
tures with the wave-vectors Kc = 0.915c* and Ka — 0.585a* 
occur in the temperature range 4.7-6 K. At temperatures be- 
low 4.7 K they vanish and an antiferromagnetic reflections 
corresponding to alternating ferromagnetic afe-planes of Ho"^'*' 



localizad moments appear. Though the spacial arrangement 
of the phases m and tT~* can not bee directly derived 
from the magnetization measurements, it is unplausible that 
the ferromagnetic in-plane interaction changes suddenly by 
switching on of the magnetic field. Therefore, we believe that 
our symbols and tt^ correspond to the real spatial se- 
quences of in-plane magnetic moments. 
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FIG. 1. Magnetic phase diagram for HoNi2B2C. H^, axis 
correspond to (110) direction. 

In this article we present a simple microscopic model for 
magnetic subsystem in the 1:2:2:1 compound which explains 
all experimental observations. We accept a model of strong 
anisotropy in which a single ion moment is directed presum- 
ably along 4 easy directions (±(1, 1,0), ±(1, 1,0) for the Ho 
and Dy compounds). Thus, the initially continuous moment 
J is reduced to a discrete variable taking only 4 values. This 
is a kind of the so-called clock model with 4 positions of the 
"hand". 

The main argument in favor of the clock model is that the 
saturation magnetization in the range of fields larger than 
7 — lOT is directed not along the field, but along the clos- 
est to the field easy direction. It means that the applied 
field us still smaller than the anisotropy field Ha- The lat- 
ter can be roughly estimated as GOT. The corresponding 
anisotropy energy for Ho^+ magnetic moment ~ 10/is is 
about 40meV ~ AOOK. Nevertheless, a single ion with in- 
teger J has no average moment in the ground state in the 
absence of external magnetic field. It happens since even a 
small tunneling amplitude w between adjacent positions of 
the "hand" (moment) leads to the ground-state in which all 4 
positions have equal probabilities. The ground-state is sepa- 
rated by a finite energy gap equal to 2juij from the first excited 
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state. A detailed analysis of the single-ion properties includ- 
ing the difference between integer and half-integer J will be 
published separately. Here we focus on the description of col- 
lective effects. 

For this purpose we introduce an angular variable 6^ at any 
lattice site r taking independently 4 values 0, 7r/2, tt, 37r/2. 
Neglecting the tunneling, the most general Hamiltonian com- 
patible with the tetragonal symmetry is: 

= i ^ [K{v - r') cos(6lr - fr) + L{r - r') cos 2(61^ - 9'^)] 

r.r' 

— hx ^ ^ cos Or- — hy ^ ^ sin Or (1) 

r r 

where K{r) and -L(r) are coupling constants and hx,y are com- 
ponents of the magnetic field. We employ the reference frame 
in which axes coincide with the easy axis directions. 

The higher harmonic terms are generated by the exchange 
interaction. Indeed, the operator of two particle permuta- 
tion specific for the exchange is a linear function of the scalar 
product of particle spins for spin 1/2. For higher spins J this 
operator contains higher powers of the same scalar product 
up to 2 J. The dipolar interaction in principle is not small 
for the Ho compound since the magnetic moment of Ho'^"'' 
is large (about 10/xs). However, it is effectively reduced in 
metamagnetic systems. The dipolar interaction between fer- 
romagnetically aligned planes is proportional to a small factor 
exp (— 27rc/a), where c is interplane and a is in-plane lattice 
constants. Returning to the exchange interaction, we find 
that for the 4-positional spins in plane only the invariants 
S1S2 = cos(((3i — ip2) and (8182)^ — (cos(^i — 952))^ are in- 
dependent, all the rest 2J — 3 invariants are equal to one of 
these two. 

The tetragonal symmetry of the single ion ground-state is 
violated spontaneously if maxisr(r) > This intuitive idea 
is supported by a variational calculation. As a result a Neel 
state arises with non-zero magnetic moments on Ho"^""" ions. 

Let us restrict the set of coupling constants to a few in- 
dependent values. We assume that the in-plane interaction 
is characterized by one nearest-neighbor negative constant K 
with all other in-plane and all in-plane Lr equal to zero. 
The in-plane interaction is assumed to be dominant to provide 
the in-plane ferromagnetic order. The interplane interaction 
is characterized by several constants Kn,L„. We shall see 
that interaction with several neighbors is essential, not only 
the nearest neighbor interaction. 

All spins in each plane are parallel. Thus, the ground state 
is determined by minimization of a spin-chain Hamiltonian: 

00 

^ [Kn cos {ifi ~ ipi+„)+Ln cos 2{(pi — ipi+„)] (2) 
z, — — 00 ,n — 1 

It should be noted that in the absence of an applied magnetic 
field it is known that the Neel antiferromagnetic state consists 
of alternating ferromagnetic a-b planes. This requirement is 
satisfied, if Ki > 0. A natural desire to simplify the model 
leaving one or two independent coupling constants unfortu- 
nately cannot be fulfilled. For example, if one leaves non-zero 
Ki and K2 and puts Li, L2 and all the rest Kn, L„ (n > 3) 
to be zero, two kind of the phase diagrams occur. The first 
diagram. Fig. 2, corresponds to < -K'2 < Ki/2 (the latter 



inequality is necessary to have the antiferromagnetic state in 
zero field). It contains 6 phases with net distribution of the 
in-plane moments symbolically depicted. Due to the symme- 
try only the sector < hy < must be considered. Fig. 3 
corresponds to K2 < 0. It is simpler and contains only 3 
phases. Neither of the phase diagrams fits the experiment 
which clearly displays 4 phases with different properties as 
shown in Fig. 1. This shows that at least the coefficients Li 
and L2 must be incorporated to describe the experimental sit- 
uation in HoNi2B2C. Moreover, we shall show later that the 
coefficient La is not zero. Thus, we restrict our model to six 
non-zero coupling constants Kn, Ln, n = 1, 2, 3. This is a 
generalization of the so-called anisotropic next-nearest neigh- 
bor Ising (ANNNI) model il- Our model differs from 
the standard ANNNI one by two features: the third-neighbor 
interaction and the 4 positions of the hand instead of two in 
the Ising model. 




FIG. 2. Magnetic phase diagram. All the parameters are 
the same as in Fig. 1 except L2 = 0.05, K3 — Lg — 0. 




FIG. 3. Magnetic phase diagram. All the parameters are 
the same as in Fig. 1 except K2 ~ —0.62, A'3 = Ls — 0. 

In order to understand why the second and third neigh- 
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bor interaction must be incorporated, one should compare 
energies of the simplest periodic sequences in the chain. The 
phases we anticipate to be realized as the ground states at 
different values of the field h are: tJ. (AF, period 2 (layers)), 
t (F, period 1), m and TT^ (period 3). Others, having 
rather close energies, are |— + (period 2); ttti, TTT^, TTi^i 
and (period 4); TTiTi, TTtn and TT^T^ (period 5); 
TTTTTi and T-IT^< * (period 6). We have found by numeri- 
cal sorting that other phases have larger energies and can be 
omitted. With these 14 phases participating in the compe- 
tition, a number of inequalities must be satisfied to ensure 
the existence of the experimentally observed phase diagram. 
Namely, on the phase boundaries the energies of the phases 
other than those being in equilibrium must be larger. For the 
reader's convenience the energies of the competing 14 phases 
are given in Table I. All are linear functions of the magnetic 
field. Therefore only their values at the corners of the phase 
diagram should be compared. 

The general investigation of the phase diagram in the 8- 
dimensional space of Kn,Ln and hx,hy is too cumbersome. 
Instead we assume that the phase diagram has 4 phase bound- 
aries, separating the experimentally established 4 phases, and 
find what are the constraints imposed by the experiment onto 
the model. The four phase boundaries found in the experi- 
ment are: 

According to the Table I, these lines are described by the 
following equations, in the same respective order as above: 

ft, = 2{Ki - 2K2 + iK:i) = i/cio; (3) 

h., + hy = 2{Ki + K2) - 4(Li + L2) = V2Hc2o; (4) 

h^-hy = 2{K^ + K2) + 4(Li + L2) = V2Hc3o; (5) 

= 2{Ki+K2). (6) 

The latter three lines intersect in the triple point /i^ = 
2{Ki + K2), hi = -4(Li + L2). Note that Eqs. (|-|) are 
equivalent to empirical equations for the transition lines found 
in Q. Thus, the theory suggests a natural explanation of all 
functional dependencies, HdiO), Hc2(0), and Hc3.{0), found 
in the experiment. The phase diagram in the plane hx,hy 
has an extremely simple shape (see Fig. 1). 

It is worthwhile to note that all above discussed functional 
dependencies were derived from purely geometrical considera- 
tions. However, the very existence of the phase diagram with 
the four phases observed in the experiment is highly nontriv- 
ial and imposes strong constraints on the coupling constants. 
These constraints are expressed as a long series of inequali- 
ties. We present here the four most important of them with 
necessary comments on their meaning: 

• Ki — 2K2 + 3K3 > 0. AF phase has minimal energy 
at hx < Hcio- 

• Ki- 2K2 + SK3 + 2(Li + 2L2)+ 6L3 < 0. The AF and 

ttl phases have lower energy than the phase tiT^< > 

on the phase boundary hx = Hdo- 

• Ki- 2K2 + 3-ft'3 + 2(Li - 2L2) + 6L3 < 0. The phases 
til and 11^ have lower energy than the phase |— > on 
the phase boundary hx + hy = y/2Hc20- 

• 2Kx - 4/C2 + 9A'3 + 4(Li + i2) + 6L3 < 0. The phases 
III and 11^ have lower energy than the phase |||^ 
on the phase boundary hx + hy — \/2Jtfc20- 



One can deduce from the second and third inequalities that: 

Ki - 2K2 + 3K3 + 2{Li +L2) + 6L3 < 0. 

Employing equations d-l) and the third inequality we obtain 
the following inequality: 



L3 < 



Hci 



Hc2l 



Hc3l 



12 



12^2 



(7) 



With the experimental values Hcio =4.1 kG, Hc20 = 8.4 kG 
and Hc3o = 6.6 kG in inequality (Q), we find the upper bound- 
ary for L3: L3 < —0.24 kG. Thus, we see that the coupling 
constant L3 cannot be zero. Thus, the experimental data im- 
ply that the interaction between magnetic planes separated 
by 3 half-periods 3c/2 is essential. Taking this interaction 
into account, we otherwise follow the principle of minimal in- 
teraction. It means that we put as many as possible coupling 
constants to be zero. In particular, we put K3 — L2 = 0. 

In the framework of our rough theory the magnetization in 
each phase does not depend on the magnetic field. It is equal 
to zero in the AF (||) phase. In the phase ||| it is directed 
along an easy axis closest to the direction of the magnetic 
field, and its absolute value is equal to 1/3 of the eeisy-axis 
saturation value. In the phase ||— > the magnetization is tilted 
by an angle arctan(l/2) = 26.6° to the easy axis closest to the 
magnetic field, and its absolute value is equal to -\/5/3 = 0.745 
of the easy-axis saturation value. In the ferro-phase | it is 
equal to 1 per site. In the experiment Q the projection of 
magnetization onto the field direction was measured. Accord- 
ing to the theory it is (1/3) cos 9h for the phase ||| (phase 2), 
0.745 cos(S — 26.6°) in the phase H— » (phase 3) and cos6 in 
the ferro-phase. While theoretical values of the magnetization 
in the phase 2 and the ferro-phase are in a good agreement 
with the experimental data, there is a discrepancy between 
the theoretical and experimental magnetization of the phase 
3 (see Fig. 5(c), in [^). In particular, in the experiment there 
is no maximum of Ms2{0) at 6 = 26.6° as the theory predicts. 
Instead the saturation magnetization decreases monotonically 
with the angle in the interval 15° < 9 < 45°. The reason can 
be that the determination of the M32 at small angle is very 
unreliable since the plateau is not clearly pronounced. On the 
other hand the values of magnetization at orientations closer 
to the easy axis where the plateau is well pronounced are in a 
good agreement with the theory. Finally, the relative differ- 
ence of the magnetization at a maximum {6 — 26.6°) and at 
9 — 45° is only 5% which may be beyond of the precision of 
the model without the tunneling taken into account {w — 0). 

From equations (^P) one can find: 



4.22 kG; K2 



1.08 kG; 



Li = -0.32 kG; L3 = -0.46 kG. 



(8) 
(9) 



Thus, we demonstrated that the low-temperature magnetiza- 
tion data and corresponding phase diagram can be naturally 
described in the frameworks of 4-position clock model withthe 
values of constants given by Eqns. (^, 

Now we consider a vicinity of the Curie temperature. We 
will show that the modulation along the c-direction naturally 
appears in the framework of the same model. The order pa- 
rameter (magnetization in a plane) is small near this tempera- 
ture allowing one to neglect the terms with the cos 2{9„ — 9'„) 
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in the Hamiltonian (|2|) (they are proportional to the fourth 
power of the order parameter s). The chain interaction Hamil- 
tonian becomes: 

oo 3 

H= ^ifnS, -S^ + n. (10) 

i— — oo n — 1 

The quadratic Hamiltonian ( |lo| ) can be represented in terms 
of Fourier-components = N^^^^ '^n-i 6"'"s„: 

H = Yk,s,s., (11) 
1 

with Kq = Ki cos q+K2 cos 2q. The value Kg has the absolute 
minimum at q = arccos(— A'i/4A'2), if \Ki\ < ^IK^l- For our 
data iCi/(4J^2) = 0.977 and q = 167° = 0.93c*. Comparing 
the theoretical value to the experimental value q = 0.915c*, 
we find the agreement to be surprisingly good, may be too 
good. We can introduce the constant to compensate a 
small discrepancy. The value K3 established in this way is 
—0.023. Though this value is not reliable, it shows that our 
minimal value was close to reality. No modulated magnetic 
phase has been found for the Dy compound [Q. From our 
point of view it means that the ratio K\/'iK2 is larger than 
1 in this compound. 

An important remark is in order: several phases which do 
not occur in the phase diagram have energies very close to the 
ground state energy. This means that a small perturbation 
(stress) can change the phase diagram drastically. 

The next step toward a more realistic theory would be to 
incorporate the non-zero tunneling amplitude w. The CEF 
spectrum numerical calculations [n4| for this amplitude give 
the magnitude w « 3 kG, which is not small, especially in 
comparison to Li and L2. The incorporation of the tunneling 
amplitude, probably weakens the strong limitations imposed 
by the inequalities. We have performed a variational calcu- 
lation of the ground-state for the extended model including 
w in the Hartree approximation. Our calculation shows that 
for zero Li the energy has the property of separability: it is 
a sum of identical functions of and hy. Therefore, the an- 
gular dependencies of Hc2 and Hc3 remains the same at least 
for small Li independently of the value of w. The variational 
calculation will be published elsewhere. 

Another important and not resolved yet question is the 
origin and the behavior of the a-modulation with the wave- 
vector 0.585a*. It appears not only in the Ho compound, but 
also in Er, Tm, and, Tb ||l^. Its wave- vector is very conserva- 
tive. Therefore, it is tempting to ascribe it to a spin-density 
wave in the conductivity electrons. This idea is supported 
by an observation of good nesting on the numerically calcu- 
lated Fermi-surface [M. However, such a treatment does not 
agree with the fact that in the Ho compound the a and c- 
modulations appear and disappear in the same temperature 
interval. Rathnayaka et al. |Q] have found an additional phase 
transition in the same temperature interval. It can be con- 
sidered as an implicit indication on the independence of these 
order parameters. From the theoretical point of view there is 
no reason for them to appear in the same point. However, di- 
rect neutron diffraction measurements do not distinguish the 
temperture where these modulations appear. 



Table I. Competing phases and their energies. 



Phase 


Energy of the phase 


AF (T L) 


— A'l + 7^2 — -ft's + il + L2 + 1/3 


F (T) 


Ki + K2 + K3 + Li + L2 + La - 




K2~ Li + L2~ L-i~ 


TU 


- ^ + K-s + Li + L2 + La - 


TT^ 


a ' a ' a a ' a 


TTTi 


L1+L2 + L3- ^ 




A'l 1 K2 1 A, aila: + lly 

2 2 2 4 


1 1 4- 


2 4 




Al 1 A2 A3 h.m+h.y 
2 2 2 4 


TTTTi 


i^ + ^ + .!i^+Li+L2 + L3-'^ 


Tun 






A'l 1 3A2 1 3A3 3Li 1 L2 1 -£-3 3/ix+2fi„ 


TTTTTi 


i(Ai +K2 + Ka) + Li+L2+La-'^ 


TiT^^ — ' 


3 + 3 + 3 3 6 
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